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FOREWORD 


This is the final report on Contract NAS8-28083. This work was in 
support of a program at the Space Sciences Laboratory, Marshall Space Flight 
Center, to study the effects of contamination on the performance of x-ray 
telescopes. 

The work performed under this contract was in two general areas* 

The first area was the modification of the x-ray ref lactometer facility to 
improve its performance and to automate data taking. This effort is discussed 
in Section 1 of this report. The second area of effort was a program to study 
the effect of contamination on x-ray telescope performance. This work is 
discussed in Section 2 and Appendix A of this report. 
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1 . 


Modification of X-Ray Reflectometer 


1. 1 Introduction 

Two modifications were made to the x-ray reflectometer 
located at the Space Sciences Laboratory, Marshall Space Flight Center. 

One was an automatic drive for the Large Micrometer Head, This system, 
interfaced with the Hewlett Packard Computer System, is used to record 
data and provides the x-ray reflectometer with an automated data -taking 
capability. Using this system, a complete scatter curve can be obtained 
automatically. Previously, it was necessary to manually reset the system 
after recording for each data point in the scatter curve. The second mod- 
ification provided an externally controlled electrical drive to move the 
microfocus x-ray source along its axis. W^ith this modification, one can 
translate the x-ray source while it is operating to locate the position 
providing maximum count rate. 

1 . 2 Automatic Drive of the La r ge Mi cr ometer Head 

As originally designed, the x-ray reflectometer sample 
detector and slit and detector assembly were positioned by a mechanical 
drive operated through a magnetically coupled vacuum feed through. After 
each measurement, the micrometer had to be manually positioned for the 
next measurement. With this modification, the micrometer is automatically 
advanced between measurements. 

Figure 1 is a sketch of the mechanical assembly. The 
motor, a vacuum rated stepping motor, drives the micrometer head through 
a six (6) to one (1) gear reduction. The motor advances 1.8 degrees per 
step. One rotation of the Large Micrometer Head corresponds to an advance 
of 25 mils. With the 6 to 1 gear reduction 1,200 steps correspond to one 
rotation of the Large Micrometer Head, Six steps correspondfj to an advance 
of . 125 mils or 1/8 mil. The unit consists of three assemblies, the Urge 
micrometer drive; the control panel; and the preset indicator. The wiring 
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schematic is shown in Figure 2. The number of steps per index is set by 
selection at the front panel of the preset indicator. In the manual mode 
(control panel switch in Disable) the unit can be run directly froizithe 
front panel of the preset indicator. Three control functions are available: 
Run (continuous stepping); Jog (advance one step); and Index (advance the 
number of steps selected on the front face of the preset indicator). In 
either the manual or automatic mode, the direction of motion (forward 
or reverse) is determined by the position of the Forward/Reverse switch 
on the control panel. 

In the automatic mode, the unit will advance during the 
computer print out operation the number of steps set on the front panel 
of the preset indicator. After completion of the predetermined number of 
steps, the control unit will be disabled for a period of time determined by 
the setting of the time delay relay in the control panel. 

This time should be selected so that the control unit is 
not enabled (after an advance) until the printing operation is completed. 

In the automatic mode, the unit will automatically advance during each 
print out until the micrometer head reaches the limit of travel in one or 
the other direction. Alternately, the advance can be limited by stopping 
the computer after a selected number of measurements. The following 
paragraphs list the instructions for installation, operation, and a listing 
of the major parts. 
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INSTALLATION PROCEDURE 


I. 

II. 

III. 


IV. 

V. 

VI. 

VII. 

vm. 

NOTE: 


Connect all cable leads to indexer unit, 

(wire markers indicate terminal strip number) 

Connect all cable leads to motor feed-thru on vacuum chamber, 
connect motor. 

Connect shielded lead to teletype: 

a. Center lead connects to lug on collector of 
QZ (power transistor). 

Top cover of the teletype must be removed 
for this connection. 

b. Shield of this lead connects to terminal #5 
of power terminal strip on rear lower left 
of teletype. 

Recheck all connections. 

Turn teletype switch to "Local” position. 

Turn on indexer unit and set controls to read 100, on dials, 

(Red "Delay" light on operators station should come on and 
after 20 second delay the red light goes out and the 
white "Waiting" light comes on. 

With teletype in Local position and inactive, adjust R, 

(lOOK sensitivity adjustment) until the red delay light 
lights, and motor indexes. Back off on R, approximately 
3/4 turn. Wait for red delay light to go out and white 
waiting light to come on. 

System is now ready to test. Depress any key on teletype, 
motor will advance, red delay light will light. After motor 
indexes through whatever number of steps you have selected 
on front panel, "Delay" then starts and will continue for 
whatever time you have adjusted on k2. 

Total time of automatic disable of motor is the combination 
of time to index plus delay time of k2 , 

This total time must be at least three (3) seconds longer than 
the time it takes to type data, otherwise an index pulse will be 
fed to the motor. 
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OPERATING INSTRUCTIONS 


I. 

II. 

III. 

IV. 

V. 

VI. 

vn. 

VIII. 


Turn on indexer. 

Wait for red "Delay” light to go out and white "Waiting” 
light to come on. 

Adjust indexer controls to number of steps desired. 

Place "Disable" switch on operator's station to "Disable". 

Type in headings and any other pre-data information. 

Pre-set detector arm by pressing "Run” button, 

("Jog" button for fine control). 

When prepared to take/record data, place "Disable" switch 
to "Enable" position. 

Motor will automatically index when first data point is typed 
on teletype. 
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PARTS LIST 


HS25V Slo-Syn Stepping Motor 
Preset Indicator, Slo-Syn SP1800-B5 
Backlash take-up board, Slo-Syn BHM144134-G1 
Brass Spur Gear, Boston Gear p/n Y2496 (1) 
Brass Pinion Gear, Boston Gear p/n YW2416 (1) 
Micro Switch 11SM3-T (2) 

Bud SC 2131 Box (1) 

MST 20 5N Switch (1) 

8-140 Terminal Strip (2) 

Potter Brirnfield Relay CDB-38-70012 (1) 

Sigma Relay 4F 8000S-SIL (1) 

IDl #1050A1 Lamp (Red) (1) 

IDI#1050A2 Lamp (White) (1) 

Arrow & Hart 81-024-GB Switch (1) 

18' Belden Cable, 6 cond. #22 (1) 

8' Belden Cable, 2 cond. #22 (1) 

Amphenol MIP Socket #77MIP8 (1) 

Ohmite Resistors #25-0362B 7, 5y\, 25 w (2) 
Ohmite Resistors #50-0565 lOOjv 50 w (1) 



1 , 3 Electrical Drive X-Ray Tube 

The x-ray reflectometer was modified to provide an 
electrical drive to translate the x-ray tube in one direction. Figure 3 and 
Figure 4 are mechanical sketches of the modified tube mount. The x-ray 
tube mount was modified to add a vacuum rated electrical motor gear box 
and limit switches and necessary mounting fixtures. This modification 
provided the capability of translating the x-ray tube along its axis using 
a control external to the vacuum system. The x-ray tube can now be 
positioned while it is at vacuum pressure and operating to maximize the 
counting rate of the x-ray reflectometer. An electrical schematic of the 
modification is shown in Figure 5. 
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MOTOR MOUNT 



"A" Holes - 10-32 Tap Thru 
"B" Holes - . 500 Ream Thru 
’'C Holes - . 187 Pin 
"D” Holes - 1/4-20 X. 70 DP 
”E" - Microswitch BA-2RU22-A2 

"F" - . 31 Diameter Hole Thru 





FIGURE 4 


MODIFICATION TO X-RAY TUBE HOLDER 




^/Z-24- T/ZA’^ffD 


SECTION 

B--B 



FIGURE 5 
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PARTS LIST 


SS25V Slo-Syn Synchronous Motor 
Bushing .75 o. d. x . 500 i, d. (4) 

Thomson External Retaining Rings C-500 
Stainless Steel Shaft , 500 dia. x 10.0” long (2) 
Micro Switch p/n BA-2RV22-A2 (2) 

Brass Gears, Boston Gear p/n Y3218 (2) 

Brass Gears, Boston Gear p/n Y3264 (2) 
Threaded shaft 12-24 
Aluminum End Plate (2) 



2 , Effect of Contamination on the Performance of X-Ray Telescopes 

Visidyne has developed a model to calculate the effect of the deposi 
of foreign materials on polished x-ray reflecting surfaces. Based on this 
model, a computer program has been developed to calculate reflectivity 
from polished surfaces and uniform thin films. The program can accom- 
odate either single layers or multiple layers of different densities or 
materials. The predictions of this model have been compared to previous 
experimental results. 

Appendix A of this document is a technical report on this work. 
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TECHNICAL REPORT 


’^Computer Code for the Calculation 
of the Glancing X-Ray Reflectivity for 
a Multilaminate Planar Surface’^ 
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INTRODUCTION 


1 . 0 

Program RFLTBL computes the reflected intensity 
as a function of the angle of incidence of an x-ray beam 
incident at grazing angles on a laminated planar surface. 

This computation is performed in two distinct steps. 

Subroutines FDELTA and FBETA calculate the real 
and imaginary components, cT and ^ , of the complex index 
of refraction from the more readily accessible physical 
parameters such as density, x-ray absorption coefficients 
and their wavelength dependence, atomic level oscillator 
strengths, and the like. Subroutine XREFL, given the thickness 
and complex index of refraction for each layer of an array, 
calculates the reflectance at a given angle of incidence. The 
main program, RFLTBL, reads in the relevant layer data, 
layer by layer, and prints it out as part of the output; calls 
FDELTA and FBETA into play to compute the required indices 
of refraction, and prints them out in tabular form; using XREFL, 
generates a stored table of reflectance values vs, angles of 
incidence and prints that out; and, if requested, passes the stored 
tables to a plotting program for direct graphical output. The 
program is a circulating one so that any number of cases 
(i, e, , a reflectance vs. angle table for a given set of layer 
data) may be processed in a single job. 
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THEORY AND APPROACH 


2. 0 


There are two major areas of computation to be 
considered — (a) calculating the reflected intensity from 
a planar laminated surface— given the necessary complex 
indices of refraction and layer depths — and (b) calculating 
the indices of refraction from more basic data. 


2. 1 X-Ray Reflection from Multi-layer Surfaces 

2. 1. 1 Electromagnetic Plane Wave in a Dissipative 
Medium 


If the charge density ^ is zero in a medium, 
and the current density is given by 


j = o“ E 

then Maxwell's equations become 


( 1 ) 


V • E = 0 


V X E 

V ‘ H 


. ^ 2. dJd 

c)t ^ at 


~ <J“E + 6 


o>E 


rel. o 


In terms of the MKS free- space electric permittivity and 

magnetic permeability, € and // , one defines the free- 

o o 

space speed of light 


c = 


7 


J-BECroWG PAGE BLANK 





and the electromagnetic field impedance 


4 = 






' Since only steady- state, monochromatic 
solutions are required, complex number methods greatly 
facilitate analysis and accordingly one defines 


E (^, t) = E (r) 


and 

H ( r, t) = H ( r ) e 

A A 

where E and H are the complex (indicated by the superposed 
circumflex) spatial amplitudes. 

In terms of these definitions, the non- zero 
Maxwell equations become 


and 


V 


V 


X 


A 

E = 


- 1 


c 






to 

c 


( e 


rel. 


- 1 




o c or 


CO 


) E 


showing that use of a complex dielectric constant 


^ = 


- i 

R I 




^ CT 


reL 


CAJ 


i* e. , where 


R ^rel. 

and 

J _ C Q- 

6a 

( and no explicit conduction current)is equivalent to treating 
a dissipative medium with a purely real dielectric constant. 
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Plane -wave solutions of the type 


i(iyt - k' r) 


and 


E (f^, t) = e 




can exist if the wave -vector k satisfies 


-A A 


k xE = - HoS) 

U> C 


and 


*k X (/h ,) = £E , 

VO ui c ^ 

A 

Inserting the first equation (for H ) into the second 

/ o tA) 

yields a condition on k : 

SO that k must be a complex quantity: 


k = kj, - i kj 

The plane-wave factor for the E- and H- fields then becomes 

A 

i(it»t - k* ^ -ki* r i{^'t - kR-"?) 

e = e e 

showing explicitly the attenuation brought about by the complex 
dielectric constant. For a pure transverse electromagnetic 
wave ( E* H = 0, kR* E = 0 = kR* H) the directions of prop- 
agation and attenuation are the same, so that 

k = (kj^ - i kj) • (kR - i kj) 

2 2 . 

- kp^ - kj - 2 i kR kj 
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which can be solved for 




-£ r ] 


In what follows it proves convenient to introduce two 
defined quantities: 


X 5 k/ {£ii/c) 


such that 


^2 = = 6 -- 


a2 
- n 


where n ( = ) is the complex index of refraction of the 
medium. Also, the non- zero Maxwell equations become 


H X E = Tl 


/\ A 


X X 'K =-£ E 
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2. 1, 2 Electromagnetic Plane Waves in a 
Multi-layer Medium 

When an electromagnetic plane wave is 

obliquely incident on a multi-layered medium, the steady- 

state result is characterized by the presence in each layer 

of two plane waves^ ^ — an advancing "transmitted" wave, 

and a returning "reflected" wave. These are, in turn, 

characterized by complex amplitudes and wave-vectors 
+ A + 

E and k— , where the + and - serve to distinguish trans- 
mitted and reflected values. At each interface of two media, 
continuity of tangential components of the E and fields 
along the interface gives rise to relationships between the 
field amplitudes and wave -vector components of the two 
media. Since each layer is treated formally in the same 
manner, a matric approach will prove most convenient. 

In the diagram below 



the pertinent variables associated with a typical n^^ layer 
are schematically indicated; where d^ is the layer thickness, 
\ is the complex index of refraction of the medium^ 
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and (p^ is the angle (relative to the interface normal) of 
the wave-vectors, k - . The unit base vector triad is defined 
as follows. The vector is the interface normal, pointing in 
the direction of the transmitted wave. The vector is a 
vector that is parallel to the interface and lies in the plane of 
incidence as defined by and k . The other vector parallel 
to the interface is defined as x , 

^h 

The two waves in the n' medium have the form 

^n 


"e- 

n 

n 


A 

i(aJt -1^ •“?) 

^ n 


i^ (ct -K“* r) 

^ ^ n 


wh« 


2 2 , 
C I 

X- r, = <a) I'n 


= n 


n 


n 


The components of K " normal and parallel to the interface 
are 


X i = X r ■ « 
"X 


— ‘ a = + n cos d) 
n X — n 'n 


and 


+ 

K - 




^ • A 

n sin © 

n I n 


For the tangential components to be continuous at all points 
on the interface, and for all times, the generalized Snell's 
Law must prevail: 

j A 

sin (p 


n. sin (b 
me I inc 


= n sin 0 ~ 

n Tn 


n+1 


n+1 


where the index of refraction of the incidence medium {usually 
air or vacuum) is taken as purely real. Since n is in general 
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complex, and n sin 6 (equal to n. sin 0 ) is real 

' n me *^inc ' 

sin may be complex: 


= “■ sin <b. / n 

* TK» r' 


n 


me 


A 

n 

me n 


in terms of which 


n cos (p n ^ n sin (p 

n n » n n ^ 


2 TT~ 

= J n n. sin Q) 

7 n me M 


me 


which also may be complex. 

In the following, two quantities occur frequently 
enough to warrant definitions 

'A ^ A 

1 I cos (p , 

_ ■ n-U = n-1 ^n-1 

n, n- 1 -C< * n cos 0 

>t n j_ n r n 


(by virtue of the generalized Snell's Law, the corresponding 
ratio of parallel components is unity, requiring no further 
definitions) 

and the phase difference, , between in-line points on 

n 

opposite interfaces ( = d ^ ) 




n cos <b 
n ' n 


To simplify the treatment of tangential field com- 
ponent continuity, it is convenient to consider an arbitrarily 
polarized plane wave as a linear combination of two plane 
waves whose electric vector is either parallel (IT-case) or 
perpendicular (jj— -case) to the plane of incidence. Contin- 
uity considerations then yield for the CT"- case; 
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= T ( 1 + a ) e 

n+1 2 n+l,n 


j ( 1 - a ) e 

n 2 n+ 1 , n 




n 


and 

e' 


n+ 1 


(1-a i (1 + a 


n+l, n 


n 


n+1, n 


) e 


-i i// 


n. 


and (in terms of H, for reasons of simplicity) for the 

iT-case: 




/A 

+ _ A r 1 x[ 

n+1 " 2 ^ ^ n 

\ n 


a ,, ] e" 

n+1 , n 


iK n 


K 

n 


1 r , / n+1 

+ 7 [ 1 - — 

2 V n 


n 


%+l,n' " 


71 


n 


and 


K 


n+1 




\2 
”n+l ^ 


1 I -| 1 ^ /-v> f 

n / n+1 , n V * ^ 


n / 
n 


+ 

n 


- i [1 - 


" \ 2 
ji+1 

A 

n 


n 


n+1 , n 


1 l-L 


n 


These multi-layer results can be expressed in terms of the 

reflection and transmission coefficients for the interface of 

( 3 ) 

two semi- infinite media : 


1 - a 


n+1, n 


n+1 ^ n 


1+a 


(q^ -case) 


n + 1 , n 


/n 

■ 


n+1 , n 


(7^-case) 


n+1 , n 


E - 

n 


n 
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2 


(cT* -case) 


n+i , n 


1 + a 


n+1 , n 


n+1 


n 


n 

jr 




( 'ft -case) 


If the equations are changed to equivalent E~ equations 

by the relationship 


2 ^ ‘ik,2^2-i2 ^ -^2 ^2 

7/ „ = (K ^ = X E - ( K • E ) = n'" 

f i n n nn nn nn 


then there obtains the simple result 


(4) 




where 



and 


M 


n 


n+1 , n 



fn 


n+1 , n 




n+1 , n 


e e’^^n 



where the appropriate lf-cd.se or €T~-case r and t 

n+l,n n+l,n 

coefficients are to be used. 


2. 1. 3 Reflection Coefficient of Multi-layer Medium 

The general geometry of the N-layer (plus 
substrate) problem consists of one semi-infinite medium (with 
real index of refraction in which both the incident beam is 

coming from a source at infinity and the reflected beam is 
returning; a set of N plane parallel layers of finite thickness 
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(with possibly complex indices n^^); and another 

semi —infinite medinm, the substrate, (with possibly complex 

index n^) in which the transmitted beam goes to infinity. The 

phase convention is adopted wherein the field amplitude in the 

k medium is defined as that at the interface of the and 
,1 ,.th 

media. The substrate, however, has no true corres- 

ponding " -1 " companion interface to its "0^”" (medium 1/ 

medium 0) interface and requires special consideration. If 

the transmitted beam, , is desired at a point 

^ i^to the substrate, a conceptual interface through d can 

o 

be considered to separate the substrate into two regions, the 
”0 " and "-1^*", both having the same index of refraction, 

n^, so tiiere is unit transmission, ihen, making the iden- 
tifications 


N+1 


E. 


me 


"N+1 

■4 

O 


E 


refl 




and using the previously derived iM-matrices: 


me 


refl 



N+1 


1m„M„ . - 

J-N iN - 1 M ^ 1 


^k \ 


"^N,k \ 
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If there is no total reflection from any interface, then the 
amplitudes Ey at the medium 1 /substrate interface are 
related to the single substrate amplitude (do = 0+) by 


so that 


K - ‘lo < 






If the matrix elements ofiM^, , are defined as 

N, 1 

Wn,i Y""" 

''“21 “ 22 / 

then the reflection coefficient, r^^, for the N -layer array is 
given by 


N E. 


“21 + “22 ■'lO 
“11 + “12 r.. 


If there is total reflection from the medium k/medium 
k-1 interface, there results 


k, k-1 k 


k-1 = 0 
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so that 


1 E, 

N, k k 


< »^N,k 


k,k-l 


r E 

k, k-l k 


^21 ■*■ *^22 \.k-l 

“ll ^-“l2'k,k-l 


N E 


Thus, in general, the reflectance R of an N-layer 
surface is given by the squared magnitude of the reflection 
coefficients, r^^, 

where is determined from the matrix elements of ftvl 

f k 

r^ k-l' k-value to be used is that of the last 

medium (in the sense N, N-l,N-2 k) to possess a non- 

vanishing reflected wave, E~. 

k 


Z. 2 Calculation of Complex Index of Refraction 


2. 2. 1 General Formulation 

Since n = i if , the calculation of n is 

essentially one of calculating £, which is defined by 




E + 4 /r p 


1 + 4 TT 
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A 

where P is the electric polarization/unit volume induced by 
the electric field E, If a sinusoidally varying electric field 
acts upon a bound electron of resonant frequency tO then 

/ C\ * 

elementary treatments show that a dipole moment^ ' 


A 

P 


e /m . 

+ i >7 ^ 

q Iq 


A 

E 


is induced, where the Vj W-term arises from radiation damping. 

If there are n such bound electrons per unit volume, then 
q 



1 + 


4lTe^ 

m 


n 

a 

- co^ + i y) ^ 

q 'q 


In separating ^ into its real and imaginary components, it is 
custornary to express the real components as the sum of a 
fixed unit term and the deviation from unity: 

i = 1 - 2 S' . 2 

so that the (reduced) real component is^^^ 


q 


2jf_l 

m 


n 


( 




(co^ 

q 


w 


) 






For frequencies far beyond the resonance 



which, as the oscillator parameters^ and Y} are absent, 

q Iq 

is os dilator -model independent. The iniaginary part. 
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since it occasions the imaginary component of the wave- vector, 
k, which in turn produces attenuation of the wave, is most 
directly associated with the linear attenuation coefficient^ 


' zV ^ 


A general treatment of € = I - 2 ^ - 2. i/^ , based 
on dispersion relations, shows that^^^ 

/to) = f J f 


<50 


/’to) = ^ 




which, following of necessity from causality, are model 
independent. Comparing the high-frequency limit of this di 
persion integral for (^) with that derived earlier 

2 


S-* 




n 


m 


60' 


Tfoi^ 




or, 


n 


m 


'rr^ e2 





This suggests defining a function (^) 


d60 

(5) 


fq(^O) 


m 


ir' 




me 


vTT 


2 2 


t 1 V.. 




to) 


such that 


«y? 




n. 
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The function ) is called the ’’oscillator strength” density 

(per unit frequency per unit volume) for the electrons associated 
with the resonance at CL} . 

q 

Accordingly, a general prescription for the real 
component, oT, is 



where the principal value of the dispersion integral has been 
implemented through the artifice of offsetting the pole at CO 
into the complex plane, analogous to the radiation damping 
term of the earlier treatment. 


2. 2.2 Specific Model for Absorption Coefficients 

Empirically it is found that the frequency 

dependence of the absorption coefficient, yCC iW), associated 

th ^ 

with the q atomic shell, is reasonably well approximated by 

( 7 ) 

the power function 

^ q'^ > c; (x ) *''’ U ( ^ ) 

q 

where 

1 th 

^ ^ shell binding energy) 


and 


Pq — empirically determined constant (of the 
order of 2 to 4) 

U (x) = 1 for X > 0 

= 0 for X •< 0 
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In terms of this model 




m c 



2 

e 


y«q(W) 


(^-^2-VPa 

17^ '^v 


and 


oo 



Two normalizing modifications to this formulation prove 
to be convenient. If N is the number of atoms per unit volume, 
then g defined as 


g = n / N 
q q 

is the oscillator strength (of the q‘^ shelD per atom, whereas 
n^ is the oscillator strength per unit volume. Secondly, 

factoring out the known 1 high frequency behavior of the 

integral yields 




?Tr ^ "NT 

^ M e •*’ e4 

m 40 2 


(p^ - 1) Re ( I ) 

q q 
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where 


q COq 


{0)^ ~ - i-n Of) 

'q 


Defining 


^ .Js. . A 3 

„ - -All 


then 


^ (generally ^ 10~^’^) 


V d V 




- i/ x^) 


Letting y = l/v 


then 


'q'Pq^*> = 2-TTT7) 


/ P„-l 


V 2 d 
1 

y ■ “T" 


±1 


Further, defining 


r = =&) 


(1 - i<) 


OC i (P - 1) 

^ 2 q 


K = 


e^ N 


( N = number of atoms/unit volume) 


J = (p - i ) I 

q q q 
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one finally obtains 


= K 


(Po J X ) ] 

H St 


where explicitly 


J (p^;x) = 


oC 

Tl i /) 


(y ) 


This integral can be integrated analytically for certain 

(rational) values of p . The early Kramers^^^ - Kallmann - 
(9) ^ 

Mark treatment assumed a value of p =3 and their results 

q 

have been used extensively. Later, examining the effects of 
neglecting radiation damping, Parratt and Hempstead^^^ 
obtained analytic integrations, correct to first order in ^ , 
for the cases p =2, 7/3, 2.5, 2.75, 3 and 4 and prepared 
tables of the results evaluated over the range 0. 01 < ^ 

q 

< 15.0. Also, they concluded that for (JO t 60 differing 

.49 

from unity by more than a part in 10 , neglect of radiation 

damping effects would produce a numerical error in the results 

^ , -3 

of about a part in 10 or less. 


2.2.3 Evaluation of the Integral J (p ;x) 

q ■ q 

The total ^ {cq ) being the sum of the partial 
(id)'s characterizing each shell, together with the sum- rule (5) 
for oscillator strengths (i . e. , Z g^ ), results in 

/to)' S aT(<») 

= k;\^ -f Re [Jq(p^;x)] 

= K 1^ l^l + Re[J^(p^;x)] - iT 


K;^ Jz+-^g Re[j(p;x) 

1 q q '■ q ' 
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2 2 
with K ■'A Z being called the normal term and K A times' 

the sum called the anomalous term. Thus, it is the quantity 

Re ( - 1 ) that is to be calculated and, neglecting the 

damping effects, this further reduces to 

I 


Jq(Pq;x) - 1 




dy 
y - ^ 


where 


T tp - 1] 


and 


r » 


- J_ - { 


to 


This real integral is singular for 0 < ^ < 1 , the singularity 
being integrable, in the principal value sense, for f 1. 

For ^ there is no singularity and a series solution is 
readily obtained: 

C<=> 




I - cxl 


h^l 


2n 


(oC+ n) 


For 0 < i, the situation is a bit more involved* For non- 


integer values of (i. e. , for not equal to an odd integer) 
series expansion of the principal value integral yields the 
result 




1 = 


OO 

t 

n=l 


1 


( - n) X 


2n 


- 'Ij' cot 7r ! 


2« 


X 


which has obvious problems when is an integer. If is 
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equal to the integer k, then letting oc approach the integer 
k as a limit gives 


J (p^;x) - 1 

<1 q 


k 


oO 

^ ^ 

_n=0 (k - n) x2n 

(n»k) 


+ 2 In x/x^^ 


Xhe programmed evaluation of the integral J - 1 proceeds 
as follows. If X is too close to unity for the series solutions to 
converge reasonably quickly, a 1000-point numerical integration 
is performed on an integrand similar to the original, complex 
one: 



_ 3 

with being set to 10 . When the integrand exceeds 100 

(i. e. , is close by the singularity), the integration step size of 
-3 , _4 

10 IS reduced to 10 . When the x- value is suitable for the 

series solutions, the series for x < I can be used without 
further consideration, but for x > 1 an additional decision 
must be made. If c< differs from an integer by more than 
one part in 100, the non-integer series is used. If aC 
is closer than that to an integer, a zeroth order approximation, 
obtained from the integer cK series, is improved with a num- 
erical o< -derivative times the difference between o< and the 
integer. The numerical o<^ -derivative is simply the difference 
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between the integral evaluated foro^ equal to the integer 
i 0. 01 (using the non-integer oC series) divided by 0. 02. 


2.2.4 Composite Materials 


2,2.4. 1 Imaginary Component — 


^^07 is the atomic absorption cross section 
for a type- i atom, N. is the number of type- i atoms per unit 
volume, and V the volume in which the atoms were found, 
then the total absorption cross-section due to all type-i 
atoms is given by^^^^ 


' O-i N, V 

- y/.v 

where = cr: is the linear absorption coefficient for type-i 

atoms. Summing over all atomic species within volume V 




mix 




cr N. V 

1^1 1 


V 


1^1 


V yU . 
^ mi 


mix 


Since y is virtually independent of the physical state of 
the material, it is a more useful quantity to deal with. 



V 


M . 
mix 



m. 

1 


M 

mix 



i mix 
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where 


m. 

f =: 

i m 


m. 


mix 


S, 

) ^ 


fh 

is the fraction (by mass) of the i atomic species. If a 
composite is known by its chemical formula — 


(ZiL _ (Z2) 


(e. g. Na^ SO^) — then in terms of the 


' ^ Hi — ^ TL'j ' Lxxcii 111 Lt^rixi 

1 Z 

elemental multiplicities, n«, and atomic weights, A : 

i i 


n. A. 
1 1 


n A. 

^ j 3 


Accordingly, for composite materials: 


mix 


4Tr 


2. 2.4.2 Real Component — 

Essentially is a measure of the polarization 
per unit volume induced in an assembly of bound charges by an 
applied E-field. For the q-shell electrons in type-i atoms, the 
expression obtained earlier for becomes 




60 

q ^ iq' OJ.^ 


) ] 


where 


K. 


2ir 2 

me 


N 


and R is the number of type-i atoms per imit volume. In 
turn, can be expressed as 

f 

N. = N P 
1 o mi; 


mix A. 
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where 


N 


f . 

mi 


mix 


Avogadro's number 


density of composite material 


fraction (by mass) of type-i atoms 


atomic mass of type— i atoms 


Thus, for composite materials 


qT • ^ c(r (Ct^) 

mix X q iq ' 


2 

e N 


2TT 




m c 


i A. 
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TESTS AND CHECK CASES 


3. 0 


3. 1 X-ray Reflection from Multilayer Surface 

The program XREFL performs both numerical 
computations and certain layer data "book-keeping" functions. 
These two aspects can be tested separately and working 
together. 


3. 1. 1 Substrate Alone 

A good overall test of the computational parts 
of the program is to apply it to the case of a substrate with no 
layers on top. The computer results obtained for = 0 

as well as several cases with r\on~zero values agreed 

exactly with tables previously prepared from the general two- 
medium interface reflectance formula^^^; 


R - f 1^ t C< - 

■ .C<- 


where 

X = e / = 9/-/^ 


and 






1 ) 


3. 1.2 Substrate and Layer Materials Identical 

A rough check of the ’’book-keeping^' aspects of 
the program is to apply it to cases where the layer material is 

pEBCroMG PAGE BI^K not FILM® 



the same as the substrate material — physically indistinguishable 
from a substrate without layers above. The results were, as to 
be expected, identical to those for the substrate alone. 

3. 1. 3 Vacuum Layer Between Layer and Substrate of 
Identical Material 

Let medium 3, the medium of the incident beam, be 
vacuum; medium 2 be of arbitrary material and thickness; medium 
1 be vacuum, but of arbitrary thickness; and medium 0, the sub- 
strate, of the same material as medium 2. Then, there is only one 
basic reflection coefficient entering in; 

n, cos <t», n cos 6 

1 ^1 0^0 


10 Ti- COS + n cos 

1 -L o o 


together with the respective layer phase shifts 

2 Jf d n cos 0 
ir = — n n n 

' n ^ 

(4) 


In terms of these, it can be shown that 


210 


^^i(i^2 - V'l) _ - yi) ^ 

^10 ~T(fTTllfV) 2 -i(^2 - V'l) -i{¥2 + ^1) -~^l) 

e - [e - e -{- e 


where 




2”5^dj 


2TTd. 




2 . 2 


- N”l ■ “l *^1 


2?rd 


27Td 




2 . 2 


sin <)>^ (by Snell's Law) 


1 


7 ^ — COS <))^ (since n^ = n^ = 1 in vacuo) 

2 It'd, 


— sin Q. 

A inc 


2 TTd 


^ — 9. (small angle approx, ) 
^ me 
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and 




2lTd, 


2 TTd 


zTTd^ 

Z'/Td 

2 

Z'fTd 


"2 ^^2 = 


2 2.2. 

n - n.- sin ©. 

' i 2 d 




2 . 2 
- sin <J) 


me 


ZlTd, 
2 


•sjn 




, 2 
sin 9. 


me 


( - I) 


A me 


2 <r. 


If the thickness d^ of the vacuum layer is taken to be 

A - % — 


2 0 . 


me 


- 1 + cos <|> 


me 


then ^ and R reduces to j *"20 j • Also, if the thickness d, 

of the material layer is taken to be 


then y =Tand R again reduces to 


d fj 0 . - 2 w _ 

me 2 


f. 


10 


I • 

Testing these results, it was found that the reflectance table 
for the layer case agreed with the table for the substrate alone only at 
those angles of incidence for which: 

2 is a multiple of'jT — and the results were the 
same for different material layer thicknesses, or, 

Z ^ multiple oijf — and the results were the 
same for different vacuum layer thicknesses. 


^ Linear Density Gradient above Substrate 

In these cases, the layer material is the same as that 
of the substrate, but of differing density. The results for 0, " 250A and 5 
50A layers are shown in Figure 1. For a wavelength of 8. 34A , 
the arbitrary (but not untypical) values of 3, 2 x lO""* and 
//oT =0.05 were assumed. In the five-layer case, the layer 

29 





densities were chosen to approximate a linear density gradient 

-^4' re respectively 0. 75, 0. 80„ 0.85, 

0. 90 and 0. 95 times the substrate density. The layer and q/' 
values, being linearly proportional to the density, varied in similar 
fashion, with the jSiS ratio, of course, remaining fixed. In the 
single “layer case, the layer density was taken as an average of the 
layer no. 5 and substrate densities, i. e. , 0. 875 times the substrate 
density. 


Qualitatively, the results are reasonable. The single-layer 
case has two pronounced interference maxima, and the five-layer 
case, a closer approximation to a linear density gradient, has but 
one maximum, substantially reduced in amplitude. In the single- 
layer case, there is a quantitative check that can be made concerning 
the separation of the interference maxima. The reflection coefficient 
for a single layer on a substrate is. given by^^^ 


21 + 


'■lO 


1 + '■21 'lO 


where and are the reflection coefficients for, respectively, 
the incidence medium/layer and layer/substrate interfaces, and ^ ^ 
is the phase shift of the layer 

hi 

y \ = ^ ”1 s 

Since the medium was assumed to have a non-zero ^ value, 
n^^^ cos (|)j will be complex. 
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where 


- z <f)^ + 4 0^ + ( 0^ - 2 ef,) 

i me 1 


2uf £- J( 9^ 

1 ’ inc 


2oC + 45/ - ( - zj": 

j* f 1 inc 1 


The phase factor in the single-layer reflection formula then b 


e comes 


-jifi 


4ltdi , 4lTd, 

""I ''r 


Neglecting the relatively slower changes in r^^ and r with respect 

10 

the maxima should correspond to Re ~ TT ^ integer* 

ztTd p 

Accordingly, for two adjacent maxima. 


oTu 

Tv R 


= 1 


a = — — 

1 2 

The single-layer maxima occur at about 0. = 94' and 0 - 

me inc 

for which the corresponding Uj^ values are 1, 374 x 10~^ and 2. 950 x lO'^ 

respectively. For A= 8.34A and Sur = 1. 576 x IQ-^, a layer 

thickness of 265A is predicted, in good agreement with the actual 

o 

value of 250A, the 6% discrepancy due probably both to errors in exact 

location of the maxima and to neglect of the angular variation on r 
and r, _ , 
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3. 1. 5 Vacuum Deposition Density Profile Studies of Parratt 

In a study of the effects of oxide layers and varying 
surface densities of vacuum deposition coatings, Parratt^^^^ presented 
some x-ray reflectance curves for copper, shown in Figure 2, 
corresponding to various multilayer approximations to certain hypoth- 
esized surface density profiles. Using the input data of his Cases I 
and III (exhibiting the most pronounced behavior) the computer program 
produced the results shown in Figures 3 and 4, agreeing, evidently, 
quite well with Parratt' s results. 

3. 1. 6 Numerical Test Point 

One final quantitative check was made by "hand 
calculating" a single reflectance value for a three-layer case whose 
parameters were chosen so as not to have any obvious "small integer" 
relationships. The hand calculation and corresponding computer 
results agreed to five significant places. 

3.2 Calculations of Complex Index of Refraction 

The programs for calculating the S and ^ components 
of the complex index of refraction likewise contain computational and 
data book-keeping aspects which can be tested separately and working 
together. 

3. 2. 1 The Integral J - 1 

q 

The central calculation in computing the complex index 

of refraction is that of the integral J (p;x) - 1. The tables of Parratt 

(7) q 

and Hempstead mentioned earlier provide an extensive and yet 
convenient check — after weeding out a number of bad entries. Upon 
encountering several perplexing, randomly occurring discrepancies 
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between table and computed values (amidst many accurate agreements), 
it became necessary to write a small program to evaluate the closed- 
form analytical expressions of Parratt and Hempstead for (p - 2, 7/3, 

2. 5, 2, 75, 3 and 4), and a table similar to theirs was produced (Table 1). 
In the few entries having four significant digits, agreement was generally 
four-place. In the two- and three-digit entries there were many small 
(last digit) and a few large (first digit) errors, due presumably, to 
extensive desk-top calculator evaluations of the rather involved 
expressions. At any rate, the series or numerical integration values 
agreed very well with the special p- value function evaluations. 

3.2,2 Values of ^ andoT* 

SincCy^^for a composite material is basically a 
mass -fraction weighted average of input data ( f values), checking 
is trivial. For oT , four -place agreement was obtained with the value 
reported by Parratt^^^^ for the pure element copper, and the value 
reported in the NAS8-27605 Final Report for the compound material 
silica glass (S|0^), thereby checking out the book-keeping aspects. 
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p = 2.0 p = 7/3 p = 5/2 p = 11/4 p = 3. 0 p = 4. 0 


1, OOOE-04 
9, 003E-04 

2. 502E-03 

1. 003E-02 

4. 055E-02 
9.286E-02 

1, 695E-01 

2, 747E-01 

4. 159E-01 

5. 039E-01 

6. 071E-01 

7. 297E-01 

8. 789E-01 
Ir 067E+00 
1.325E+00 

1. 740E+00 

2. 252E+00 
2. 354E+00 

1,950E+00 
1.674E+00 
1. 531E+00 

lo439E+00 
1, 373E+00 
lo 324E+00 

1« 254E+00 
1.207E+00 
1, 147E+00 

1. 098E+00 
1, 059E+00 
1, 039E+00 

1, 013E+00 
1„ 003E+00 
1. OOlE+00 


2.405E-03 

9. 471E-03 

1. 727E-02 

3.608E-02 
6. 061E-02 

5. 857E-02 

2. 226E-02 
-5. 681E-02 
-1. 909E-01 

-2. 845E-01 
-4.010E-01 
-5. 464E-01 

-7. 310E-01 
-9. 728E-01 
-1.311E+00 

-1. 868E+00 
-2. 559E+00 
-2. 722E+00 

-2. 195E+00 
-1. 840E+00 
-1. 657E+00 

-1. 541E+00 
-1.459E+00 
-1. 397E+00 

-1. 310E+00 
-1.252E+00 
-1. 178E+00 

-1. 119E+00 
-1. 071E+00 
-1. 047E+00 

-1. 016E+00 
-1. 004E+00 
-1. OOlE+00 


2. 056E-03 
9. 543E-03 
1, 884E-02 

4. 445E-02 

8. 976E-02 
i. 120E-01 

9. 918E-02 

3. 921E-02 
-8. 364E-02 

-1. 760E-01 
-2. 950E-01 
-4, 481E-01 

-6.467E-01 
-9. 121E-01 
-1. 289E+00 

-1. 917E+00 
-2.698E+00 
-2. 895E+00 

-2. 308E+00 
-1. 916E+00 
-1.714E+00 

-1. 587E+00 
-1. 497E+00 
-1. 429E+00 

-1. 335E+00 
-1. 272E+00 
-i. 192E+00 

-1. 128E+00 
-1. 076E+00 
-1. 051E+00 

-1. 017E+00 
-1. 004E+00 
-1. OOlE+00 


1. 399E-03 
8. 051E-03 

1. 758E-02 

4. 794E-02 
1. 156E-01 
1.704E-01 

1. 934E-01 
1. 666E-01 

6. 802E-02 

-1. 797E-02 
-1. 362E-01 
-2. 959E-01 

-5. 114E-01 
-8= 081E-01 
-1. 240E+00 

-1. 972E+00 
-Z, 889E+00 
-3. 144E+00 

-2. 469E+00 
-2, 021E+00 
-1. 794E+00 

-1.650E+00 
-1. 549E+00 
-1. 474E+00 

-1. 369E+00 
-1. 299E+00 
-1. 210E+00 

-1. 140E+00 
-1. 083E+00 
-1. 055E+00 

-1. 019E+00 
-1. 004E+00 
-1. 002E+00 


9. 210E-04 
6.311E-03 

1.497E-02 

4. 595E-02 
1.271E-01 
2.082E-01 

2. 653E-01 
2. 747E-01 
2.071E-01 

1. 320E-01 
1. 960E-02 
-1.414E-01 

-3.682E-01 
-6. 914E-01 
-1. 174E+00 

-2. 008E+00 
-3. 062E+00 
-3. 380E+00 

-2. 619E+00 
-2. 119E+00 
-1. 866E+00 

-1. 707E+00 
-1. 596E+00 
-1. 514E+00 

-1. 399E+00 
-1. 323E+00 
-1. 227E+00 

-1. 150E+00 
-1. 089E+00 
-1. 060E+00 

-1. 020E+00 
-1, 005E+00 
-1. 002E+00 


3. OOOE-0-^ 
2. 698E-03 

7. 481E-03 

2. 970E-02 
1. 151E-01 
2.449E-01 

3. 987E-01 
5.440E-01 
6. 308E-01 

6. 288E-01 
5.775E-01 

4. 561E-01 

2.325E-01 
-1. 468E-01 
-7. 897E-01 

-2. 004E+00 
-3. 606E+00 
-4. 225E+00 

-3. 141E+00 
-2. 448E+00 
-2. 107E+00 

-1. 895E+00 
-1. 750E+00 
-1. 643E+00 

-1, 495E+00 
-1. 398E+00 
-1.278E+00 

-1. 183E+00 
-1. 108E+00 
-1.072E+00 

-1. 024E+00 
-1, 006E+00 
-1. 002E+00 


Table 1 
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PROGRAM OUTPUTS 


4. 0 

The program always returns printed output, and optionally 
returns plotted graphical output. 

4. 1 Printed Output 

The printed output consists of two main sections. In 
the first section, the input data is reproduced with appropriate 
labelling. In the second section, the calculated results are presented 
in tabular form, 

4.1.1 Input data 

The contents of the case label card are printed at the 
top of the page. Just below will be the general run parameters; 

1. the in vacuo incident wavelength 

2. the number of layers above the substrate 
(always layer no. 0) 

3. the number of incidence angle / reflectance 
value tabular entries to be calculated 

4. the minimum incidence angle to be used 
in the table 

5. the incidence angle increment per tabular 
entry, and 

6. the choice of plotted return or not 

Then, layer by layer (topmost layer with highest index, substrate 
with zero index), the layer parameters are presented: 

1. the depth or thickness of the layer 

2. the layer material density, and 

3. the number of atomic elements of which the 
layer material consists. 
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For a given layer, the parameters of each of the constituent atoms are 


successively presented; 

K the atomic number 

2, the atomic weight 

3. the fractional (by mass) abundance of the element 

in the layer material 

4* the mass absorption coefficient, f > ^.t the 
specified incident wavelength 

5, the number of atomic shells sufficing to describe 
the inner atom, and 

6, the atomic shell parameters 

a) the shell index ( NB a single index can refer 
to one level, e. g* , Lj, or to a group, 

e* g* , JVIt 

b) the shell absorption edge wavelength 

c) the exponent of the power -law dependence 
of the shell's contribution to the total 

f , and 

d) the oscillator strength (per atom) of the 
shells 


1. 2 Calculated Results 


The contents of the case label card are printed at the 
top of the page, and just below, in vacuo incident wavelength is repeated. 
This is followed by a table of the layer parameters required for calcula- 
ting the reflectance: 


1) the layer thickness 

2} the computed value, and 


3 ) 


the 


:omputed^^- 


value. 
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This, in turn, is followed by the table of reflectance values consisting 
of 

1) the angle of incidence ( in arc minutes) 

2) the (same) angle of incidence (in radians) 

3) the reflectance value, and 

4) the stopping layer 

The stopping layer is defined as that layer in which the total 
attenuation (due both to inter-layer reflection and intra-layer absorption) 
of the incident beam field strength exceeds 1 Total reflection by a 
layer is interpreted as the beam stopping within that layer. At zero- 
incidence angle, there is no penetration (similar to the case of total 
reflection) and the top-layer index is arbitrarily assigned as the 
stopping layer. Finally, in the case of the substrate, the value -1 
is used to signify net transfer of beam power into the substrate (whether 
subsequently absorbed or transmitted) and the value 0 to signify total 
reflection from the substrate. 

4. 2 Plotted Graphical Output \- 

When the plot index is set to I (a plot is desired), the 
incidence angle /reflectance value table is transmitted to a plotting 
program package (drawing upon standard CALCOMP software) consisting 
of GOPLT, which is a plotter initializing routine; FRMPLT, which 
determines internal scaling values, draws a graph frame with tick 
marks and labels it; CRVPLT, which plots the reflectance vs. angle 
table; and ENDPLT which closes out the plotting routines. 

Since each computation center has its own plotter conventions, 
GOPLT makes a call to a subroutine PLTID3, which transmits two 
arguments — a six-character BCD word (to identify user's plotted out- 
put) and an overall scale factor (a value of 1 , 0 producing output with 
an 8 1/2" x 11" format as indicated by the bounding tick marks). The 
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subroutine PLTID3 is best written by st systems programmer who 
knows the computation center conventions and should perform the 
following tasks: 

lo make a call to PLOTS to define the output device 
and a PLOT working buffer 

2. space the plotter away from previous user*s plot 
results 

3. label the plotter paper with the present user-s 
six- character BCD ID word 

4. define a new origin: 0* 5'^ or more away from 

user label and 0. 5" up from paper^s edge, and 

5. RETURN to GOPLT. 
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5. 0 INPUT DATA FORMAT AND DECK SETUP 

The general data deck consists of the data for an arbitrary 
mber of cases (each case being a reflectance vs, angle curve for 
g' en set of layer material data) followed by a ''=i-*=;-END OF DATA"^!''!'-''' 
The structure of a general data deck is given in Figure 5, The 
structure of a partial deck for a single case is as follows; 


I. Case Label Card 


Cols. 

F ormat 

Item 

1-72 

12A6 

Arbitrary case identification 
data (e. g. , date, specimen 
description, etc) 

General 

Cas0 Paratnctcr 

Cols. 

F ormat 

Item 

1 - 10 

ElO. 3 

In vacuo incident wavelength 

14 - 15 

12 

(Angstroms) 

Number of layers above substrate 

18 - 20 

13 

( 0 or more - up to 20) 
Number of reflectance table 
entries to be calculated ( 1 or 

21 - 30 

ElO. 3 

more - up to 300) 

Initial incidence angle value 
(arc min. ) 

31-40 

ElO. 3 

Incidence angle increment 
(arc min. ) ' 

45 

11 

Plotted result index ( 1 = YES, 


0 = NO) 

followed by layer description subdecks ( 1 or more - up 
to 20), containing 
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***END OF DATA*** 


Case #2 

i 


- 

Layer //3 



(Case if 2 Arb, Info) 

(etc. for additional cases) 


IqLAM (1,5) 8 PQ(1,5) i GQ(1,5) I QLAM(1,6) i etc 
|qLAM(l,l) i PQ(1,1) 8 GQ(1,1) i QLAM(1,2)I etc, etc 
Z (l) T AuTI FM(1) I FMURHO(l) I NQ(1) I 

^e.g. 7 shells) 

D(3)l RH^i I NATMSl 

^(e.g. = 1: layer #3 also pure element) 


■r 

Layer #2 


QLAM(1,5)I PQ(1,5) I GQ(1,5)| 
iQLAMa.l) B PQ(l,l)i GQ(1,1) I QLAM(1,2) i PQ(1,2) 
Z(l) S~A(1) S FM(1) B FMURH0{2) B NQ(1)1 

^-<e.g. 5 shells) 


D(2) I RH0 I NATMS 

i(e.g. = 1: layer Ifl is pure element) 


etc, etc 



f f |qlAM( 3,5) i PQ(3,5) 

/ |QLAM(3,l) i PQ(3, 

/ "”|Z(3) i A(3) g FM(3) i FMURH0{3) § NQ(3) 

/ ^(e.g. 6 shells in atom #3) 


j Atom if 2 QLAM(2,1) § PQ(2,1) “ GQ(2,1) 

' 2(2) @ A<2) i FM(2) FMURH0(2) iS NQ(2) 

c #1 ^(e.g. 2 shell in atom #2) 

f |qlAM( 1,5) fi PQ(1,5) B GQ(1,5) 

Atom #1 |qLAM(1,1) £1 PQ(1,1) S GQ(1,1) R QLAM(1,2)B PQ(1,2) 


(e.g. 3 atomic species on layer #1) 
(e.g. CuSO. ) 


ACLAM i NLAYRS i NPTS 1 THMIN B DTHINC If IPL^iT I 
^(e.g. 3 layers) 

(Case #1 Arb. Information) 


FIGURE 5 










III. General Layer Parameter Card 


Cols. 

Format 

Item 

1 - 10 

ElO. 3 

Layer thickness (Angstroms) 

11 - 20 

ElO. 3 

Layer material density (g/cm ) 

24 - 25 

12 

Number of chemical elements in 



layer material (lor more — up 



to 10) 

followed by element de 

scription subdecks ( 1 or more — 

up to 10), 

containing 


Atomic Element Parameter Card 

Cols. 

F ormat 

Item 

1 - 5 

F5. 2 

Atomic number of element 

6-10 

F5. 2 

Atomic weight of element (g/mole) 

11 - 20 

ElO. 3 

Element abundance (mass fraction) 

21 - 30 

ElO. 3 

Mass absorption coefficient 



^ (cm /g) 

35 

1 1 

Number of atomic shells ( 1 or 


more — up to 8) 

followed by 


V. Atomic Shell Parameter Cards ( 1 or 2) 


Cols 

• 

F ormat 

1 - 

8 

E8. 2 

9 - 

13 

F5. 2 

14 - 

18 

F5. 2 

19 - 

26 

E8. 2 


Item 

Shell 1 absorption edge (Angstroms) 
Shell 1 absorption coefficient power 
law exponent (positive) 

Shell 1 oscillator strength (per atom) 
Shell 2 absorption edge (Angstroms) 
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27 - 31 


F5. 2 


Shell 2 absorption coefficient 
power law exponent (positive) 

32 - 36 F5. 2 Shell 2 oscillator strength 

(per atom) 

37 - 44 E8. 2 Shell 3 absorption edge (Angstroms) 

etc. , etc, 

NB No second atomic shell parameter card is 
required if the layer material consists of 4 or less 
atomic elements. 

In order that computer system run termination 
diagnostics not be printed immediately below the last 
output value table (but rather on a new page), the last 
card in the data deck must be one containing the 
characters ***END OF DATA*'!'* In colums 1-17. 
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6.0 PROGRAM LISTING 

The program listing is given in the following data 

sheets. 
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o n o o o 


PRO(iHAM RPLTdU (iNPUTtOUTPUT) 

MAIivI program to COMPUre the glancing x-kay 
angle Of" iNClOtNCL FOR A LAMlNATpo SURFACE* 

for”each [AreS! kekraction 


HEFLECTIVITr WS« 

the beta and delta 

ARE CALCULATED 


dimension title a^) 9NLR(2D oO( 20) 9^(10) ,A(10> .FMIIOl.FMiiBUftMn* 

pRim 300 
son FORMAT (IH I ) 

IHLT s 0 

I read loOt (TITLEU) tl = 102) 

100 rOHMATdSAoJ 

IE ( title U ) sNE« GH'iMtttENO) GO TO 5 
IE ( IPLT«NEe 0) CALL ENUPLT 
EXIT 

101 EOKMAT(eiO«JOX«T292X9i3,2ElO,3»4X.n ) 

Te-?^r»*^S l^HPLOTl ( IsyES^OSNO) S 9 I 7 ) 

TF(IPLOT.NE,0) iHlT s iplT ♦ I ^ 

N[Rr='NU4;«Tl‘'’'-“'’‘^“-'’ 60PUT(AHI0 »RD,1.0) 

DO 20 IL = i»NLRl 
NLR(IL) s! NLRl - II 
read 102» U ( IL ) 9 RHO 9 NATMS 

102 EORMAT (2Ei0.3,3X» 12) 

20 a 22 ex. 7 H 0 .PTH . 

• G2CH3.2eX, ’ 

DO 10 lA s 1»NATMS 

.03 F0RHA;?2%f2!?^;S”t:;:y;r’’^’^"'’”“”‘*'^“” 

• 2?0 3H L’"FTr2%”':i«":lt^ P*P*M.TFRS.,’.::/“;{Si.3HZ ..FS.O. 

• RH.lImOA) -^fIo P«ACT. ='.E10.3210X.6HmU2RH0. 

NSHELS = JidM CM22S.2lOX.UHNO. AT. SHELLS ..12) 

lOA F0RHAT(4(E»“2f2F5!2n’’'’°“*’'®’’®“*'*’ISI.IS.l.NSHELS) 

2^4 FoilJiTut^x;i7H?TSMK*jiRu"p;L*«?^^®“'"’J5U 
20 
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PRINT 205» (TITLEU) 

205 FORMAT (1H1«5X,12A6) 

IF (NUAYRS,NE.0> go to 30 
PRINT 206» YACLAM.NLAYRS 

206 F0HMAT(//6X»13HLAMDA(VAC.) »iE10,3*10H ANGSTROMSf //6X»5HREFL, » 

♦ 18H parameters of TH£*I3*IIM tAYERS. . • , »//0X*5HLAYER*3Xt 

♦ 9MTHlCKNESS,7X.5HOELTA,9X,4MBETAt/9X*3HNO.f3X»llH( ANGSTROMS) , 

♦ //20(9X,I2*3E1A,3*X)> 

PRINT 207 

207 F0RMAT(7X*6H(N0NE> »//lX> 

GO TO 40 

30 PRINT 206. VacLAM.NLAYRS. (NLR(I) .0(1) .OELTA(I) .BETA(I) ♦ 

♦ IsI.NLAYRS) 

4C PRINT 208. OELTA (NLRI ) .BETA (NLRl ) 

208 format (6X.21HAN0 THE SUBSTRATE. ••,.// 1 OX . IHO 1 14X. 2E14. 3> 

PRINT 209 

209 FORMAT (/////6X.31HTA8LE OF REFLECTANCE VALUES. ////12X* 

nHTHETA(lNc.).lOX.UMREFLECTANC£.3X.8HSTOPPlNG./6X»5H(ARC . 

♦ SHMIN.) .3X.9H(RA0IANS) »2lX»5MLAYER./lX) 

NPTSI r NPTS ♦ I 

DO 50 N a 1. NPTSI 

TrilNCM(N) a ThMIN ♦ FLOAT(N • 1>«DTHINC 
THINC * 2.90888£-4«THINCM(N) 

CALL XREFL<THINC.VACLAM.NLAYRS.0. delta. beta. R2(N) .NREFL) 

50 print 210. THINCM(N) »ThINC»R2(N) .NREFl 

210 F0RMAT(6X.F5.1.2E15*3.I9) 

PRINT 200 

iFdPLOT.EU.O) GO TO 1 

THMAX a THMIN * FLO AT ( NPTS ) *DTHlNC 

Call frmplt (Thmin.thmax.o.-o.i.i.i.thincm.rb.nptsd 

CALL CRVPLT(THINCM,R2.NPTSi,IH »n) 

GO TO 1 
END 
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FUNCTION FUELTA<VACLAM.RHOiNATMS.Z,A»FM»NO,OLAMfPQ,GQ) 

^UO) »A(10) ,FMU0l eNOdO) sOLAMdOjB) 9 PQ (10*8) *60(10*BJ 

MoUM • UoO 

no 20 lA = 1*NATMS 
QSUm s 0.0 
NO I A » NO I I A) 

00 iO lO = 1»NQIA 
IF(G0(IA*IQ) .EO.O.n) go to 10 

10 CONTlNUt^*^"’ * (P(J(IA*IO) *0LAW{IA*IQ)/VACLAM) 

20 ASUM S ASUM ♦ FM ( I A ) ♦ ( Z ( I A ) ♦ 0SUM)/A(IA) 

FOELTA - 2*70200E-6*HHO*VACLAH*<»2#aSUM 

return 

FNC 


FUNCTION F8£TA(\/ACLAM,RH0*NATmS.FM*FMURH 0) 
UiMENSlON FM(IO) *FMURHO(10) 

SUM s 0.0 
00 10 I B 1*NATmS 
10 SUM = SUM ♦ Fm{ I) »FMUHHO( I ) 

FBETA s 7.95775E-10 *VACLAM*RhO*Sum 

return 

end 


DO 

10 

?0 


FUNCTION FJULl (P*X) 

IF(P,gT«K0) go to 20 
P R i i\i J ^ o 0 

FORMAT (////23H *«**« p i$ LESS ThAN 
FJOLl s O.U 


It////) 


return 

(X.GT.1.CE3) GO TO lO 
FJOlI c -l.o 
IF (X.lT. 1 0 OE-?) return 
ALPH B (p - l,0)/2,0 
/I a 1.0/X*«2 

IF(X.LT.O.S) FJOLl = FJXLO(ALPHtXl) 

IF (X.GTtltl) FJQ*-1 s FJXhI ( ALPH t XI ) 

FJQU = FJNUM(ALPH* 

ENU 


XI) 
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FUNCTION FJNUM(A#Xn 

F(Y> a (Y - AI)»Y*»A/({Y - Xl)*»2 ♦ l.OE-6) 

YN » 0.0 
OY = l.OE-3 
PINT = 0.0 
DO 20 ISUm « 1*500 
FNO = F(YN) 

FNl s F(YN ♦ OY) 

FN2 a f<YN ♦ 2.0«DY) 

IF(AMAX1 {ABS(FN0> ?ABS(FN1) »A8S(FN2) ) .GT.1.0E2) GO TO 30 
Fl»^T a PINT ♦ FNO ♦ 4.0*PN1 ♦ pN? 

10 YN * float (ISUM)*2.0*rJY 
20 CONTINUE 

FJNUM = A*FINT*DY/3.0 - 1.0 

return 

30 FINTl a 0.0 
YNl a YN 

UO 40 iSUNl a 1*10 

FNO a FUND 

FNl a FCYNI * 0.1*DY) 

FN2 s F(YNl ♦ 0.2«>DY) 

FINTl a FINTI ♦ FNO ♦ 4.0*FNl ♦ FN2 
40 YM s YN ♦ float USUMD <»0.2*OY 
FINT a pint ♦ FINT1*0.1 
GO TO 10 
END 


FUNCTION FJXL0(A*XI) 

SUN a 1,0/ ( (A + 10.0) *XI) 

DO 10 I a 1,9 

10 SUM a (SUM ♦ I.O/IA ♦ FLOATUO - I)))/XI 
FJXLO a -A*S>Um - 1.0 

return 

END 


FUNCTION FJAHI (AtXI) 

FN a IFIX(A + 0,b) 

EPS - ^ aa pfV 

IF (ABs(EPS) .LT.0*01) go TO 10 
FJXHI a FJA(A,XI) 

return 

10 FJXHI a FJN<FN,XI) ♦ (FJA(FN ♦ 0,01, XI) - FJA(FN - 0.01, XI)) 
♦ *£PS/0.02 

RETURN 

ENU 
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function FUN(FN*XI) 

FJN * -1,0 

IF<FN,gE,11,0) return 
SUM » XI/(FN - 26.0) 

DO 10 I « 1|25 
FNLI b FN - FL0AT(26 • n 
rF<A8S{FNLl) .6T.0.01) 60 TO iO 
SUM s XI»SUM 
IF(I.E0.2&) go to 20 
U' SUM « XI*(SUM * 1.0/FNLI) 

20 continue 

FJN s FN«(SUM - XI»*FN*AL06<XI) ) 

return 


function FJA(A.XI) 

DfiTA Pi/3, Ui5927/ 

SUM a Xi/(A - 26.0) 

DO 10 I s 1»2S 

10 SUM s X!*(SUm * 1,0/ (A - FLOAT (26 - I))) 
ARO s PI«(AHO()(A»1,0) - 0.5) 

FJA * A*(iUM . XI*»A*PI*^IN(AR6)/C0 S<aR6) ) 

return 

fno 
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SUBROUTINE *HEFL ( THJNC* VACLAMfNU YR »0 , DELTA »BETA* RMA62 » NREFU 
OiMENSlON U<20)» l>ELTA(21)* BE T A < 2i) , cm ( 2 , 2) * CML(2*2) 
COMPLEX CM* CML* FNCSL* FNCSLl* FNCOS* RL» TL* OL* R 

if(thinc.ne.o,o) go to 10 

RMAG2 *1.0 
NREFL * NLAYR 

return 

10 S (1. 0*0.0) 

CM(l*2) s (C. 0,0.0) 

CM(2,1) s (0. 0*0*0) 

CM(2,2) s (1. 0*0.0) 

NLRl * NLATR ♦ I 

D(NlR1) s 0.0 

TNET s 1.0 

DO 30 NL a 1 *NLR1 

NREFL * NLRI - NL 

FNCSLl s CmPLX (TMInC.0.0) 

if<nl.gt,i) Call sfncos<thinc*delta(nl-d *beta(nl-i) ,fncsld 
Call sfncos (Thing * delta (nl) *beta (nd *fncsd 

RL a (FNCSLl - FNCSU/(FNCSL1 ♦ FNCSL) 

TL a (2.0*0.0)»FNCSU/<FNCSLI ♦ FNCSL) 

OL a CmPLA(6,283185*0(NL)/VaCLAM, 0.0)*FNCSL 
IF(CA8S(TL) .LT.l*0E-4) 60 TO 40 
IF (-AIMA6 (DL) .gT*9,0> 60 TO 4q 
TNET * TNET*CABsUl>*EXP(AImA6(Di ) ) 

IF(TNET.LT.1.0E-4) 60 TO 4o 
IF (NL.NE.NLRl ) 60 TO 2o 

IF(REAL(FNCSL)/CAeS(FNCSL) .6T.l.nE-4) NREFL * -1 

60 TO 30 

20 CML(1*1) a CEXP ( (0.0* 1 .0) *DL) /TL 
CML(1*2) a NL«CEXP( (0.0,«1.0)*DL) /TL 
CML(2*l) * RL»CEXP( (0.0*1. 0)*DL)/TL 
CML(2*2) a CEXP( <0.0*-l.O)<^DL)/TL 
CALL MaTmPY (CM*CML*CM) 

30 continue 

40 R a (CM(2*1) ♦ RL*CM(2*2) )/(CM(l,l) ♦ HL*CM(1,2)) 

RMAG2 a CABb(R)*»2 

RETURN 

ENU 
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subroutine SFNC0S(THINC*DELTA>8ETAiFNC0S) 
complex fncos 

FTUb = SQWT (0,5* (SORT ( (THINC«*2 > 2,0*DELTA) * 4 *0<»BETA**2) 
« ^ ABS(TMINC**2 - 2.0*DELTA))) 

TF(9£TA,NE«0.D) Go to 10 
FNCOS a CMPLX(FTDB*0«0) 

JF(THINC**2.LT,2»0*OELTA> FNCOS a CMPLX (0» 0 * -FTOB) 

return 

10 FNUOS = CMHtX(FTUB,-ABS{BETA)/FTnB) 

IF (THINC**2#LT.?.0*DELTA) FNCOS a CMPLX ( ABS (BETA) /FTOS*-FTOB) 

return 

FNU 


SUBROliTlNE MATMPY (Cl tC2*Cl2) 

OlMEiMSiON Ci(2,2)» C2(2*2)» C12(?*2). C(2»?) 

Complex ci» c?, Ci2, c 

DO 10 I=l«2 
DU 10 J=l»2 
CU»J) S (0.0, 0,0) 

DO 10 K*1»2 

lU CU 9 J) s CU*J) * Cl ( UK) »C2(K, J) 

DO 20 JSU2 
DO 20 J=1.2 
20 C!2(UJ) a C(UJ) 

return 

END 
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SUHrtOUTiME GOPLT (PRI0»FACT) 

XTOT = 510.0 
YTOT = 11.0 

CALL PLilD3(PRID»ATOT,YToT,FACT) 
call plot (‘J. ?5».b»3> 

CALL PLOK'^.TS* . ^»P) 
call plot (S».7?^.7.H,3) 

CALL PLOT(^.75»8.0f?) 

CALL PLOT (0.<95» H*U»2) 
call PLOT C^.S5»H. 0.2) 
call PLOT{0.0».b»-3) 

return 

FiNU 


SUBROUTINE CRVPLT ( X * Y ,NtBCU» IC IRC ) 

COMMON/PLTR/XMAXL,XMlNC»YMAXCfYMlNC»DX»DY»XINCHf YINCH 
DIMENSION X(Jno) tYOOO) 
dimension AMOoOOO) »YMOt)<300) 

IC s 0 

IFdCIRC .NE. 0) IC s I 
II = 0 

DO I B 1* N 

IF(XMINC ,LE. X(I) .AND. X<I) ♦Le# XMAXOII » II ♦ 1 
IF (XMINC »LE. X(l) .ANU. X(I) .lE. XMAXC) XMODdD ■ Xd) 
TF (XMINC .LE. Xd> .AND. Xd) .lE. XMAXC) YMOD(II) * Yd) 

10 continue 

IFdl .EQ. 0) RETURN 
DO 20 I B 1, II 

IF (YMOUd) *LT. YMINC) YMOQ d ) * YMINC 
IF (YMOUd) .GT. YMAXC) YmOD(I> b YMAXC 
20 continue 

CALL I INE (XMOD* YMODf 1 1 ♦ I . IC* 1 ♦XMTNC»DX»YMINC*DY» .ofl) 

XWHKRE = (XMODdI)-XM!Nt:)»XlNCH/ (X maXc-XmINC) ♦ •! 

YWHERE - (YM0D(U)-YMINC)*YINCH/(YMAXC-YMINC) -.05 
CALL symbol (XWHERE fYWHEREf .l»BCOtO.O»l ) 

return 

FNO 
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SUBROUTINE FRMPLT (XMlNt XMAX tNYCALC* VMIN» yMAX ♦ X » Y » N) 

comnon/pltn/xmaxu»xminc*ymaxc*yminc.ox*dy,xinch,vinch 

OrMENSlUN X (300) »Y (300) 

DATA XINCH/9,0/*YINCH/7*0/ 

CALL PLOTdl *0»0i&*-3) 

IFIMYCALC *E(J, 0) UO TO 20 
YMAX s Y(1) 

YMjh) s V (1) 

DO iC I * 1» N 

IF (x(i) ,lt. xmin .or. x(I) .gt. xmax) go to 10 

TF (V(I) ,Lt, YMlN) YMIN s Y(l) 

IF (Y(I) ,GL. YMAX) YMAX » Y(I) 

io continue 

20 IF (XmaX ,LE. XmIN) return 
S aBS(XMAX - XN1N)/10.C 

IF(ALOGlO(OX) ,LT. 0.0) NMiNXsIFiX (ALOGlO (DX) ) -1 
IF (ALOGIO (UX) ,L». 0.0) GO TO 3() 

NMI^x § ALOGIo(DA) 

30 CONI INIJE 

nXRU S ux*l 0 .0«'<» <«=NMINX) 


IF 

(DXRU 

eGE . 

l.U) 

OX 

a 

1.0 

If 

(OXRO 

.GE. 

2.0) 

OX 

s 

2.0 

IF 

(OXRO 

.GE. 

2.5) 

DX 

s 

2.5 

IF 

(OXWU 

• G£ . 

5.0) 

OX 

s 

5.0 

OX 

a OX* 

10.0»«NMINX 





XWNFRM = float UFIX (XMIN*1 .OOOl/OX) )*OX 
If (XMIN ,L1. 0.0) XMNFRM s XMNpRM - QX 
XMXFRM s float ( IFIX (XMAX#0.9999/DX) )*DX 
IF (XMAX ,GI. 0.0) XMXFRM = AMApfiM * OX 
IF (YMAX ,Lt. YMlN) RETURN 
DY s ARSIYMAX c YMIN)/io,0 

IF ( ALOGl 0 (OY) .LT. 0.0) NMInY s IFIX ( ALOGl 0 (OV) ) ®1 
IF (ALOGlO (DY) ,Ll. 0.0) GO TO 
NMJNY = ALOGIO(UY) 

40 CONTINUE 

OVNO s OY»10.0** (-NMINY) 

IF (DYRO «(jE. 1.0) OY s i,o 
IF(0YRD .G£. ?,o) DY = 2,0 
IF(DYPn sGE. 2,5) DY = 2.5 
IF (DYPO .6E. b.O) OY a 5.0 
UY s DY*10.0 **NMINy 

YMNFRM = float ( IF I X ( YMIN«l .000 l/pY) ) «OY 

IF (YMIN ,LI. O.O) YMNFrm s YMNFRM « DY 

YMXFRM a float { IFJ X (YmaX*0 .9999/oY) ) «OY 

IF ( Ym^X ,Gr. O.O) YmXFRM a YmXfPM ♦ dy 

XMINC s XMNFrm 

KMAAC s XMXFRm 

Vising a ymnfpm 

YMAXC = YMXFHm 

OYX s (XMAXC • XM I NC ) I 0 . 0 / ( OX*X inch ) 

OVY s (YMAXC - YMlNCd*10.0/ U)Y«»YTNCH) 
nx s (XMAXC - XMINC)/XINCH 
DY = (YMAXC - YMINO/YINCH 
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CALL AXIS (O.OtO.O*2lHlNC. AN6l£ (ARC MIN, ) » -21 #XlNCH*O.Of XMlNCf 
♦ OXtUVX) 

call Axis (0,0*0. 0*11HREFLECTANCE,U*YINCH, 90. 0 *YmINC*DY*0VY) 
CALL PLOT( 0,0*YINCH*3) 
call PLOT(XINcH,YInCH*2) 

CALLi PLOT(XINCH,0,0*2) 

IF(XMaXC*XMINC ,<3E, 0.0) GO TO flO 
XZERO B -XMlNc^XlNcH/ (XMAXC-XMIN c) 

CALL PLOT (XZERO*0,0#3) 

CALL PLOT (XZeRO*YINCM*2) 

80 continue 

IF (YMAXC*YMInC .GE» 0,0) GO TO 90 
YZERO B -itmINC*YIncH/(YmaXC-YmINC) 

CALL PLOT ( 0,0*YZEKO*3) 

CALL PLOT (XINCH,Y2EP0,2) 

90 continue 

XTEAR a XINCH ♦ .?S 
XTEARl = XINCH ♦ ,95 
XTEAR? a XINCH ♦ ,55 
YTEAR s YINCH ♦ 

YTEAHI s YINCH ♦ ,5 
CALL PL0T(XTEAR, 0.0*3) 

CALL PL0T(XrEAR*,2,2) 
call PLOT(XtEAH* Y teAH,3) 

CALL PLOT(XTEAft* YTEAttl *2) 

CALL PLOT (XrEARI*YTEAHl*2) 

CALL PLOT (XrEAR2*YTEAm*2) 

return 

ENu 
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